LPTHE 00-18 
May 2000 



On the Consistency of the Exact 
Renormalization Group Approach Applied to 
Gauge Theories in Algebraic Non-Covariant 

Gauges 



M. Simionato 

LPTHE, Universite Pierre et Marie Curie (Paris VI) et Denis Diderot (Paris 
VII), Tour 16, l er etage, 4, Place Jussieu, 75252 Paris, Cedex 05, France and 
Istituto Nazionale di Fisica Nucleare (IN FN) Frascati, Italy 
E-mail: micheles@lpthe.jussieu.fr 



Abstract 

We study a class of Wilsonian formulations of non-Abelian gauge the- 
ories in algebraic non-covariant gauges where the Wilsonian infrared 
cutoff A is inserted as a mass term for the propagating fields. In this 
way the Ward-Takahashi identities are preserved to all scales. Never- 
theless BRST-invariance in broken and the theory is gauge-dependent 
and unphysical at A ^ 0. Then we discuss the infrared limit A — > 0. 
We show that the singularities of the axial gauge choice are avoided 
in planar gauge and light-cone gauge. In addition the issue of on- 
shell divergences is addressed in some explicit example. Finally the 
rectangular Wilson loop of size 2L x IT is evaluated at lowest order 
in perturbation theory and a noncommutativity between the limits 
A — ► and T — ► oo is pointed out. 
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1 Introduction 



In recent years, a lot of activity has been devoted to the study of Quantum 
Field Theory in the framework of the Exact Renormalization Group Equa- 
tion (ERGE) of Wilson [|TJ . The principal advantage of this approach is the 
fact that the ERGE can be solved numerically at the non-perturbative level: 
this feature in principle pave the way to a number of important applications. 
However, in practice, the method works well for scalar theories, where there 
is an extensive literature and a collection of reliable numerical results (see 
|2| for a recent review and a complete bibliography) whereas its application 
to gauge theories should be regarded as problematic. The root of the diffi- 
culties becomes clear if we remind the logic of the ERGE approach in the 
formalism of the one-particle irriducible (1PI) effective action |J, which 
is the following: i) take a given Quantum Field Theory and introduce a fic- 
titious infrared cutoff A by modifing the bare propagator in some way; ii) 
write the functional identity which describes how the 1PI effective action of 
the theory changes when the infrared cutoff is lowered; iii) fix the boundary 
conditions at the ultraviolet scale Ao by identifying the effective action with 
the local bare action and integrate down numerically || or analytically in 
perturbation theory M the equation up to the the infrared value A Ao- 

This logic is perfectly legitimate in scalar theories, but quite delicate in 
the case of gauge theories since there is no way of inserting the infrared cutoff 
in a way consistent with gauge-invariance (we mean with BRST-invariance): 
as a consequence the A ^ theory is gauge- dependent and unphysical and 
the limit A — > must be taken strictly. This generates two distinct problems. 

The first problem, which has been extensively studied in the literature 
[|5], [(J, is the following. Since the infrared cutoff breaks BRST-invariance, 
the usual Slavnov- Taylor identities have to be modified and assume a very 
complicate form || such that it is practically impossible to find truncations of 
the effective action consistent with the modified identities. As a consequence 
at the non-perturbative level there is no rigorous control on the fact that the 
gauge-invariant theory is recovered when the infrared cutoff is removed. On 
the contrary at the perturbative level this is guaranteed by a theorem: no 
matter how the infrared cutoff is inserted, the A — > theory is gauge invariant 
provided that we choose the ultraviolet bare action in a very precise way, 
which can be established in perturbation theory to all orders by solving the 
fine-tuning equations discussed in ||. Therefore, at least in principle we could 
consider solved this problem, even if in practice the fine-tuning equations are 
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quite involved and cannot be solved out of perturbation theory. 

The second problem, which in our knowledge has been treated only par- 
tially in the literature, is the question of the regularity of the A — > limit. 
Actually the smoothness of this limit is guaranteed only for the proper ver- 
tices at non-exceptional configuration of momenta f7]: but these are not 
physical quantities. On the contrary, the limit is extremely delicate in the 
computation of physical quantities as for instance the interquark potential: 
we will see in this paper that the limit must be taken with great care in 
order to have sensible results. In such a situation, a naive numerical analysis 
of the problem is certainly questionable. Apparently, this subtle point has 
never been recognized in the previous literature and is explicitly analyzed 
here for the first time. 

The guideline of our research program consists in setting up a formalism 
where these problems can be faced analytically with perturbative computa- 
tions, in such a way to have definite answers on the recovering of the gauge 
invariance and on the regularity or the singularity of the A — > limit for 
physical quantities. In this logic, in a couple of recent papers || |9| we pre- 
sented a Wilsonian formulation of gauge theories in which the Wilsonian 
infrared cutoff A is introduced as a mass-like term for the propagating fields. 
In this way the Wilson's Renormalization Group Equation for the 1PI ef- 
fective action becomes consistent with the Ward-Takahashi identities to all 
scales. This important property, which allows to solve the problem of finding 
gauge-consistent truncations but not the question on the regularity of the 
A — > limit, was proved in the Abelian case in covariant gauges || and 
in the non- Abelian case in algebraic non-covariant gauges ||. In this latter 
case, the crucial point was the fact that in algebraic non-covariant gauges the 
ghost fields decouple and the gauge-symmetry can be implemented via sim- 
ple Ward-Takahashi identities instead that in terms of complicated Slavnov- 
Taylor identities. However, as we said in ||, even if the implementation of 
the gauge symmetry is quite efficient in these formulations there are various 
disadvantages to be taken in account. The first disadvantage, common also 
to the covariant version of the method, is the need for an explicit regulariza- 
tion of the evolution equation, since when the mass cutoff is employed the 
ultraviolet momenta are not sufficiently suppressed. However, in our view 
this is a little price to pay since, at least in perturbation theory, there exist 
gauge-invariant techniques to manage the ultraviolet problem (dimensional 
regularization, higher derivative regularization, etc.). The second disadvan- 
tage is the gauge- dependence of the would-be physical quantities at A ^ 0. 
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However, even this problem is not so dramatic since the gauge-dependence 
can be controled throught the generalized functional identities described in 
||. It is only the third disadvantage, i.e. the fact that non-covariant gauges 
have a very delicate A — > limit, which is a very serious one. This feature 
gives rise to difficulties since if the A — > limit does not exist for quan- 
tities which are instead well defined in other gauges the whole approach 
should be considered as suspicious. In practice this fact does not prevent 
phenomenological applications since there are infrared safe quantities which 
can be computed in some approximation (i.e. by a truncation of the evolu- 
tion equation) without encountering problems. Nevertheless in questions of 
principle one would like to have a formalism with a well definite perturbative 
expansion and fully consistent with other gauges for any physical quantity. 
This is why in this paper we carefully study the infrared limit of the non- 
Abelian theory in different algebraic non-covariant gauges. We will see in 
particular that the planar gauge and the light-cone gauge are regular in the 
A — > limit whereas the axial gauge is singular. The origin of the problems 
of the axial gauge will be connected with the spurious infrared divergences 
introduced by the double pole term in the propagator, i.e. the term propor- 
tional to n 2 q^q v /[n ■ q} 2 . This is hardly unexpected, since it is well known 
that in the standard approach using the Cauchy Principal Value (CPV) pre- 
scription the singularities connected with the double pole are responsible for 
many serious problems, as the loss of unitarity [10| and the failure of the 



Wilson loop consistency check |TT|. In particular this latter quantity is the 



most interesting to study since it is connected with the interquark potential, 
which is the typical quantity one would like to compute non-perturbatively 
in the Wilson renormalization group approach [|T^]. We recall that problems 



are expected to be present in the perturbative computation of this quantity 
since it was shown long time ago, working in the temporal axial gauge fI3 



that the Wilson loop of sides 2L x 2T which directly corresponds to the in- 
terquark potential computed at order g 4 is different from the corresponding 
computation in covariant gauges. More dramatically, in the standard per- 
turbative expansion, the temporal gauge Wilson loop does not exponentiate 
in the limit T — > oo i.e. it is not of the form 

H/ rLT T =°°expH2T\/(2L)], (1) 



sec 



n[] for a critical analysis). The same problem was found in the spatial 
axial gauge for loops in the Euclidean space with the side T in the direction 
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of the axial vector [] [fTIJ] . The origin of the problem was found in the 
Cauchy Principal Value Prescription which in fact breaks gauge- invariance. 
There has been a lot of work in the literature to try to solve this problem 
(see the references cited in the monograph ||15|| ). In the Wilsonian approach 
one naturally obtains a regularization of spurious divergences different from 
the CPV prescription and which is consistent with unitarity. Nevertheless, as 
we anticipated in |J, still one expects to have problems in the infrared limit 
A — > 0. These problems are discussed in detail in this paper: in addition 
we prove that they are absent in more stable algebraic non-covariant gauges, 
as the planar gauge and the light-cone gauge. This is not a surprise since 
it is well known in the literature that these gauge are safe with traditional 
prescriptions (the CPV for the planar gauge and the Mandelstam-Leibbrandt 
(ML) for the light-cone gauge) and in particular the Wilson loop consistency 
check is successfully passed [TR, [H| at least at order 0(g A ) in perturbation 
theory. The light-cone gauge in particular should be considered as a promis- 
ing starting point for our analysis. However we stress that to prove that 
the Wilsonian i.e. massive formulation tends to the standard massless for- 
mulation (which is safe) in the A — > limit is a subtle point since there 
are infrared sensitive quantities where the massive version and the massless 
version of the same gauge choice are different, i.e. the A — > limit is singu- 
lar. The non-trivial point is to show that this does not happen for physical 
quantities. In principle this could also be false in a finite order perturbative 
computation. As we said, one of the motivation of this paper it to set up a 
workable formalism where such a question can be answered with a reasonable 
analytic effort. 

In order to compare and contrast with the previous Wilsonian literature, 
we remind that a Wilsonian approach to algebraic non-covariant gauges has 
already been presented in [T^, [T7| ; nevertheless our analysis differs from that 
approach, where a generic cutoff function is employed thus totally obscuring 
a number of important properties and difficulties of the formalism: in par- 
ticular the gauge-invariance issue and the infrared problems can be clearly 
studied only within our framework. In spite of the peculiarities of our ap- 
proach, we stress that our results should be regarded as quite general, since 
there is a large class of cutoff functions which in the infrared behave as a mass 
term: for instance the most used cutoff in the numerical ERGE literature || 

1 Obviously there exist loops which give the same result in axial gauge and in Feynman 
gauge and in general smooth loops are expected to be safe. 
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the exponential cutoff 



R A{q) = f^fif'l (2) 
1 — exp(— q 2 /A 2 ) 

is in this class since R\(q) — A 2 for q 2 <C A 2 and thus soft momenta are 
effectively screened as if there were a mass A 2 . This is not, however, the 
general case, since for instance the sharp cutoff is not in this classQ. In 
general we expect that the situation of the infrared limit for other cutoff 
functions is worse than in the mass-cutoff case, since the gauge-symmetry 
is completely broken, Ward-Takahashi identities have a cumbersome form, 
and possible gauge-cancellations between different divergent contributions 
are lost. Actually one could speculate that at the non-perturbative level, due 
to infinite resummations of classes of Feynman diagram the infrared problem 
could be eventually solved. However, as a matter of fact, the perturbative 
expansion is the only controled analytical tool we have at our disposition. 
This is why we are interested in giving a Wilsonian formulation in which i) 
a well defined perturbative expansion can be defined and ii) it is possible 
to perform practical, analytical computations with an effort not too much 
bigger than in the standard massless approach. In particular in this work we 
discuss a general method and a series of tricks to analytically compute finite 
parts of one- loop Feynman diagrams in the A — > limit. This is relevant in 
preparation to a computation of the Wilson loop up to order 0(g A ). 

The plan of the paper is the following: in section 2 we review the basics 
of the standard gauge-fixing procedure, pointing out some subtle aspects 
which are usually not noticed since an infrared regularization consistent with 
BRST symmetry is assumed, which is not our case. In section 3 we study the 
properties of the general linear gauge L fl (p)A'^(p) = in presence of a mass- 
like cutoff. In particular the massive axial gauge, the massive planar gauge 
and the massive light-cone gauge are analyzed. In section 4 we discuss the 
problems of the axial gauge as connected to the presence of the double pole 
in the propagator, whereas in sections 5 and 6 we show that these problems 

2 This can be seen for instance by using the approximation to the sharp cutoff 

q 2 cxp(-(< Z 2 /A 2 ) fo ) 



Ra(q) 



l-exp(-(g 2 /A 2 ) b ) 



with b — > oo. fn this case the infrared modes are very strongly suppressed since R\(q) is 
divergent as q 2 — > 0. The same is obtained with other regularizations of the sharp cutoff. 
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are absent in planar and light-cone gauge. In section 7 the Wilson loop for 
a finite rectangular path is explicitly computed up to order 0(g 2 ) and it is 
shown that the physical limit A — > exists and is independent of the gauge- 
fixing choice. Section 8 resume our work and contains our conclusions. Two 
appendices close the paper: appendix A is a simple but instructive exercice 
that shows the gauge-dependence problem in the computation of a physical 
quantity, the pressure of a free photon gas; appendix B collects the technical 
tools and tricks which are used in the text to compute one-loop Feynman 
diagrams. 

2 Remarks on the gauge-fixing procedure 

In a strictly gauge-invariant theory one cannot define a gauge-field propaga- 
tor since strict gauge-invariance implies that the free two-point function is 
transverse and therefore cannot be inverted. As a consequence the starting 
point of the perturbative analysis in terms of Feynman diagrams is miss- 
ing. Thus, in order to define a perturbative quantum field theory from a 
classical gauge invariant field theory one is forced to break gauge invariance 
through the addition of a gauge-fixing term to the classical actionP[ The way 
gauge-invariance is broken depends on the chosen gauge-fixing term, which 
is at large extent arbitrary: the essential requirement is that the gauge-fixing 
term, depending on the gauge field and other fields (ghosts, antighosts and 
auxiliary fields), must be a BRST-cocycle: in this hypothesis it is possible to 
prove that the physical quantities are gauge-fixing independent and therefore 
that there are not physical arbitrarities. 

In the path integral approach the gauge-fixing is usually introduced by 
means of the Faddeev-Popov argument: since there is an infinite number 
of gauge-equivalent configurations giving the same contribution to the func- 
tional integral the partition function is ill defined. Then the problem is solved 
by integrating only on representative of any gauge orbit. The representative 
are fixed by means of a gauge fixing condition of kind jF a (A) = 0; typical 

3 Outside perturbation theory, for instance in the lattice approach, one can define the 
theory without an explicit gauge-fixing. However in order to show that the continuum limit 
exists, which has been rigorously proved only in perturbation theory a gauge-fixing term 
is needed even in the lattice formulation. In general the gauge-fixing term is unavoidable 
to make contact with the perturbative formulation, which is the only one we analytically 
control. 
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choices are JF a (A) = d^At (covariant gauges) or J ra (A) = n^A^ (algebraic 
non-covariant gauges). These latter are interesting even because they avoid 
the problem of Gribov copies, i.e. the choice of the representative of the 
gauge orbit is unique |0| . 

Now it is convenient to introduce a bit of notations (the full list of our 
notations is reported in 0). A gauge transformations of parameters 5uj = 
5u ■ T = 8uj a T a is denoted by 

5 G A M = D^Su, 5 G ip = ig5utp, 5 G i) = -igtp5uj, (3) 

where T a are the generators of SU(N C ) in the fundamental representation 
and N c = 3 is the number of colors. The set of physical fields of theory is 
denoted by = (A"(p), i[>i(p), and the bare action with >Sb(0); the 

expression |£ denotes the differential operator jj-D^ . A gauge-invariant 
ultraviolet regularization, which exists in non-chiral theories, is understood. 
With these notations the gauge-fixed partition function can be written as 

Z T = j [dcf>] det 5 -£ 5(F(A)y s ^\ (4) 

The usefulness of the gauge-fixed approach is based in the following funda- 
mental property: if the observable 0(0) is gauge-invariant, then its mean 
value in the gauge ^(A) = 0, defined as 

< 0(0) > T = ±- J [rf0] det 6 -£ 5(^(A))e^^O(0), (5) 

is independent of the gauge- fixing function J ra : 

5 G O(0) = O <0>r=<0>j», T^T'. (6) 

This means that in the computation of physical quantities we have the free- 
dom of choosing any possible gauge-fixing. 

We think that it is useful to review here the formal argument that support 
this statement with the scope of showing where and why this argument fails 
in the Wilsonian formalism. 

The first step consists in introducing anticommuting fields C and C and 
commuting auxiliary fields A in order to rewrite the determinant det ^f- as a 
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functional integral over Grassmannian variables^ 



det oc / IdCdC] exp 
Ouj J 



-iC ■ 5 -fc 

OUJ 



(7) 



and the functional delta function as 



8{F{A)) = J [dX] exp [iX ■ J 7 (A)] ; (8) 

then the partition function (defined up to a multiplicative constant) can be 
rewritten in the form 

Z T = J [dCdCdXd<P] exp [iSjr(<f>, C, C, A)] (9) 
where the total bare action 

- dT 

Sf = S B + X- HA) - C ■ g^D^C (10) 

is invariant under the BRST-tranformation 

sA^ = DyP, sC = -^gC xC, sC = A, sX = (11) 

for any choice of the gauge-fixing function. The cohomological property of the 
BRST tranformation s 2 = allows us to prove the independence of BRST- 
invariant quantities from the gauge-fixing choice. Consider for instance a 
class of gauge-fixing functions ^{A) determined by one or more continuous 
parameters e: then we have to prove that physical quantities are independent 
of e. This can be easily shown by observing that the gauge- fixing term can 
be rewritten trivial cocicle 

A ■ T e {A) ~ C ■ = s{C- T £ {A)) ■ (12) 

Using the fact that trivial cocycles do not contribute to the partition function, 
i.e. 

J [rf$]e^ ($) s /($) = 0, (13) 



4 Here we use the hermiticity conditions C = and C = — [Ell such as the relation 
(C ■ f£C)t = C ■ C holds. 
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(we have denoted by $ = (0, C, C, A) the full set of fields of the theory) one 
sees that the partition function is e— independent 



^Z £ = J [d$]e iS **h (c ■ ^e{A) \ = (14) 

and the same holds for BRST-invariant observables 0(<£>) such as sO = 0: 

^ < O > e = y J [d^eWh (o{*) C ■ ^s(A) \ = 0. (15) 

Unfortunately this proof is purely formal: a rigorous analysis must face the 
question of infrared divergences which can invalidate the argument. In fact, 
the perturbative expansion has no meaning at all, if we do not specify care- 
fully as the infrared divergences are managed. The need for an infrared reg- 
ularization can conflict with BRST-invariance and in general to show that 
gauge-invariance is recovered when the infrared regulator is removed is non- 
trivial since the existence of the infrared limit for physical quantities is a very 
delicate point. To our knowledge this issue has been not enough emphasized 
in the literature on the Wilson Renormalization Group, in spite of the fact 
that it is a very crucial one. Actually, even if the recovering of the gauge 
symmetry at A = has been proved for proper vertices at non-exceptional 
configurations of momenta j7|, there are no theorems guaranteeing the in- 
frared safety of important physical quantities as the interquark potential, 
which is the first thing one would like to compute in the Wilson renormal- 



ization group approach [[12]. In particular, as we said in the introduction, 
there are problems in its perturbative computation. The reason is that the 
fundamental property of the exponentiation of the Wilson loop at T — ► oo 
is lost at A ^ since it is related to highly non-trivial gauge-cancellations 
between different contributions which are lost when BRST-invariance is bro- 
ken by the infrared cutoff. In other terms we expect a non-commutativity of 
limits T — > oo and A — > 0. This can be seen even at 0(g 2 ) in perturbation 
theory using the non-covariant formalism and it will be discussed at length 
in section 7. This peculiarity suggests that numerical computations of the 
interquark potential should be performed with great care. 
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3 The general linear gauge 



In general we are interested in linear gauges where J 7 {A) is a linear function 
of the gauge field, 

F{A) = Z/(p)A», L» = arv + &p„. (16) 

Notice that this is a class of interpolating gauges involving both the covariant 
gauges (if a = 0) and the algebraic non-covariant gauges (if 6 = 0). Other 
non-covariant gauges, as for instance the Coulomb one, can be obtained for 
particular values of a and b. In order to integrate out the auxiliary fields it is 
convenient to introduce a small convergence factor in the functional integral 
measure by adding to the BRST action a term 

i&A • A, 

where £2 is a positive parameter of mass dimension —2. This terms is a 
trivial cocycle and therefore does not change the essential properties of the 
system; in particular by using standard techniques it is possibile to prove 
that physical quantities are formally independent of £2- By using the linear 
equations of motion 



A a = -^rL^A 11 (17) 
£2 



we can eliminate the auxiliary fields A a and consider the reduced action 

Sun. g . = J ~\f,u -F^-C- L»D»C + j i^A, ■ L v A v . (18) 

To this reduced three-level action, we add an infrared cutoff as a mass-like 
term 

SJStgXA C, C, A) = S& B XA, C,C) + J \a 2 A\ (19) 

Now we can define an invertible massive propagator for the gauge fields. The 
explicit form of the propagator can be obtained in Euclidean space (Euclidean 
notations are recalled in appendix B) by inverting the matrix 



(20) 
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where t^ u {j> E ) denotes the transverse projector in Euclidean space 

VW=V-^?. (21) 
Pe 

In particular we are interested in the £2^0 limit, when one eigenvalue of 
the propagator vanishes, due to the transversality property 

L^D A ^(p E ) 6 =° 0, V A (22) 
and therefore D\^ v is not invertible (for any A). The final result is 



n ( \ 6 - 1 x _ P ■ l ( l i*Pe,v + L u p E ,y) 
Va ^ [Pe) ~ p 2 E + A 2 ^ (p| + A 2 )((p E ■ L) 2 + L 2 A 2 ) + 

Pe^Pe,v L 2 A 2 L il L v 

(p| + A 2 )((p E ■ Lf + L 2 A 2 ) (p 2 E + A 2 )((p E ■ L) 2 + L 2 A 2 ) ' 



(23) 



Notice that in limit £2 - > the propagator is invariant up to rescaling of 
the gauge-fixing L{p E ) — > C{p 2 E )L{p E ). In particular for oc n M we recover 
the massive axial gauge introduced in || which satisfies simple Ward iden- 
tities, whereas for oc p^ we obtain a massive version of the usual Landau 
gauge, which does not satisfies linearly broken Ward identities, but instead 
satisfies a little modification of the usual Slavnov- Taylor identities (this is 
the Curci- Ferrari model in Landau gauge [0). We have also computed the 
explicit form of the propagator at £2 7^ 0, but the resulting expression is not 
particularly illuminating and there is no scope in writing it here. 



3.1 Planar gauge 

The standard massless version of the planar gauge is formally obtained from 
the generalized axial gauge if we replace the parameter £2 with the momentum 
dependent function ^(p) = n 2 /p 2 ; in this case the BRST action reads 

Si*Z = J ~\f,u -F^-C- n»D,C + A ■ n% + ^A • ^A (24) 
or, after elimination of the auxiliary fields, 

S r p tar = J ~\f,u -F^-C- n»D,C - inM„ • d 2 n»A u . (25) 
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The planar gauge is less problematic than the axial gauge, but the price to 
pay is in a more complicate realization of symmetries: in opposition to the 
axial gauge case, the action (p5| ) does not satisfies simple Ward-Takahashi 
identities (the reason being the presence of derivatives in the gauge-fixing 
term) but instead Slavnov- Taylor-like identities; moreover ghost fields play 
a non-trivial role even if less crucial than in covariant gauges^. The planar 
gauge has been studied quite intensively in the literature and has some inter- 
est in itself, and also in comparison with the axial gauge and the light-cone 
gauge in the problem of inconsistencies of perturbation theory. Moreover 
the computation of Feynman diagrams in planar gauge is similar but simpler 
than in light-cone gauge. For these reasons we will analyze in detail this 
gauge choice. We found that the more convenient way to insert the infrared 
cutoff in order to have a simple propagator is to modify the massless BRST 
action as 



S p lanar(A) = J —F^ ■ F^ v - C ■ D ,C + A • nM„ 
This gives as tree level reduced action 

s£L(A) = J -\f,u ■ F*» + h 2 A ■ A + 

- ^n^" • (d 2 + A 2 )n^ - C ■ n^D^C . 



(26) 



(27) 



Notice that the "mass" A 2 multiplies the term 

\a» [g^ - ^) A u = \x 9l] A\ i,j e {0, 1, 2} 



and thus only transverse (with respect to n^) degrees of freedom are screened. 
There are other possible ways of introducing the infrared cutoff in the planar 
gauge, but this is the more interesting one in the sense that one obtains a 



5 For instance the planar gauge has the interesting feature that all diagrams involving 
ghost loops identically vanishes ^5|. This property is preserved even in the massive 
version. Nevertheless, there are non- vanishing vertex corrections including ghosts as ex- 
ternal lines. In principle this could be avoided by considering a formulation without ghost 
fields p3[ but then the Ward identities becomes quite cumbersome. 
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propagator which is simple and very similar to the propagator of the light- 
cone gauge. Its explicit form in Minkowsky space is 



-D A Jp) = 1 L v _ ± n »P» + A2n » n » | . (28) 
^ v p 2 — A 2 + is \ ^ [p ■ n] [p ■ n] 2 J 

The most important difference with respect to the massive axial gauge case 
is the fact that in the massive planar gauge there are spurious divergences at 
p-n = even at A ^ 0. Therefore we need an explicit prescription to manage 
them. The simplest choice which works, at least up to the order 0{g i ) in the 
stardard massless casef], is the CPV prescription 

1 , p ■ n , „ 

- lim p - . (29) 



[p ■ n\ e^o (p ■ n) 2 + £ 

We will use this prescription even in the massive case. We stress that if we 
split the propagator in two pieces as 

Da,hv = Da,iiv + D\^ u (30) 

where D\^ u is the term proportional to n^n^/n 2 , 

° A ' M = -[p.n] 2 (p 2 -A 2 + tey (31) 

we can see that our formulation reduces to the usual planar gauge expression 
in the A — > limit if the Da^u term can be neglected. This is obvious in the 
computation of infrared finite Feynman diagrams, far from being obvious in 
the computation of infrared divergent Feynman diagrams, and definitely non- 
trivial for physical quantities sensitive to the infrared, such as the interquark 
potential. 



3.2 Light-cone gauge 

The light-cone gauge is the most used and the most tested algebraic non- 
covariant gauge. At the present it passed many serious consistency checks 
and should be considered at the same level of safety of the covariant gauge, 
at least for what concerns the perturbative expansion in four dimensional 

6 Actually in higher order computations the CPV prescription could be problematic; 
see the discussion in section 5.1. 
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gauge theories. In particular it has been explicitly proved at order 0(g 4 ) in 
perturbation theory that the Wilson loop computed in the light-cone gauge 
gives the same result of the Wilson loop in covariant gauge, provided that we 
use the Mandelstam-Leibbrandt (ML) prescription to regularize the spurious 
singularities fl9 l. 

The Euclidean massive light-cone propagator can be formally extracted 
from espression ( ^3|) for light-like vectors = Ue,^ such as n 2 E = 0: 

D '^ = ^+Tf\ s - ^ feT^h (32) 

The light-cone condition can be realized in Euclidean space if the gauge 
vector tie = (ni, n 2 , n 3 , 714) has the form n E = (0, 0, 1, i) or n* E = (0, 0, 1, — i), 
or an equivalent one after a 5*0(4) rotation. Writing 

p 3 =p||cos0||, p 4 =P\\ sin 6\\, p"\\ = Pe ■ n E p E ■ n* E = pi + p 2 4 , (33) 

we obtain the explicit expressions 

1 1 -ipA+P3 PE-n* E e~ ie w 



Pe ■ n E iPa +P3 Pl+ p\ p\ P\\ 



(34) 



After Wick rotation in Minkowsky space po = —ip4 we see that this approach 
corresponds to regularize the spurious poles with the ML prescription. The 
Minkowskian propagator reads 

_n n- 1 / ^uPv + ^vPp, , A 2 n fl n u \ 

Uk^W - p2 _ A2 + ie ^ [[p . n]] + ii p . n]] 2j> ^ 



with 

1 _ p ■ n* 



[\p ■ n]] (p ■ n){p ■ n*) + ie 



n= (1,0, 0,1), n* = (-1,0,0,1) (36) 



and reduces to the standard one at A — > if 'the D\^ u term can be neglected. 
Differently from the planar gauge, the light-cone choice is extremely conve- 
nient since not only it avoids the double pole problem but still mantains 
the advantages of the axial gauge, i.e. the transversality property which 
garantees the full decoupling of ghost fields, 

n»D^ A (q) = 0. (37) 



14 



Furthermore the simple Ward identities 

q " D ^ [q) = ii^j] (38) 

and 

ifDf (q)V^ p (q, p, -g - p)Df (g + p) = (g) - D* T (g + p) (39) 
hold, where 

^ P (p, g, r) = (g - r) M ^ p + (p - g) p ^ + (r - p)^ (40) 

is the non-Abelian three-gluon vertex. By using these ingredients we explic- 
itly checked the transversality of the gluon propagator and in general the 
Ward identities on the proper vertices with direct diagrammatic considera- 
tions, for any value of A. We repeat that in planar gauge, on the contrary, 
the Ward identities are much more complicated and the gluon propagator is 
not transverse with respect to or p M . 



It is important to stress that: i) it is impossible to impose condition fl39|) 
in covariant gauges: this is the reason why the gauge symmetry is unavoid- 
ably more complicate in the covariant case and ghost fields have to be taken 
into account ii) differently from planar gauge, we checked that there are no 
other possibile forms of the light-cone propagator consistent with (|37|),(|38|) 
and (|3"9"|): the transversality constraint plus the symmetry requirement fixes 
univocally the form (|35| ) for the light-cone propagator. In other words, equa- 
tion (|35|) is the unique way of introducing a Wilsonian infrared cutoff in a 
non-Abelian gauge theory consistently with Ward identities. 



4 Problems of the pure axial gauge 

As we recalled in the introduction, the standard approach to axial gauge 
with the CPV prescription has various problems and in particular the test of 
exponentiation for a properly chosen Wilson loop fails. There has been a lot 
of work in the literature to solve this problem (see for instance the citation 
list in [jTjJ ) , but in spite of this effort in our opinion at present there are no 



completely satisfactory solutions. For instance there is an apparently simple 
solution consisting in changing the prescriptions on the gluon propagator and 
the ghost interaction in such a way that the Wilson loop becomes the same 
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as computed in covariant gauges. This is the logic of Cheng and Tsai KB 



and also of the approches based on interpolating gauges, in which one try to 
define the axial gauge as a deformation of a more regular gauge. This kind 
of approaches has a long hystory (starting from the old work of jl3) until 
very recent papers [^7], |28fl ) nevertheless they are not completely satisfactory 
since i) the transversality property n^D^ = is lost; ii) the ghost fields are 
no longer decoupled; iii) the Ward identities have no more a simple form; iv) 
a careful study of the infrared singularities appearing in the limit in which 
the modified gauge tends to the pure axial gauge is needed. For such reasons 
these modifications are so drastic that in fact they should be interpreted 



as a switch to a truly different gauge [15|. The situation in the Wilsonian 



approach is different but still there are problems which will be discussed in 
the following sections. Our final conclusion will be that the axial gauge is 
definitively sick, at least in perturbation theory, and alternative gauge choices 
should be considered. 



4.1 The CPV regularization of the double pole 

In the Cauchy Principal Value prescription, the double pole is defined as the 
derivative of the single pole, 

± = -±± (41) 

where l/[ps] is defined as in (^). This means that the double pole is regu- 
larized as 

1 p 2 3 - A 2 

= &(p| + A2) 2 • (42) 

We stress two things: 

1. One could expect some problem with this regularization of the double 
pole, since we loose the property of the Euclidean propagator of being 
positive definite. This can be seen by looking at the eigenvalues of the 
Euclidean pure axial gauge propagator (at zero mass and with the CPV 
prescription), which are 
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The first three eigenvalues are obviously non negative; the problem is 
with the last one, which is not positive definite. This means that the 
integral 

w) = L w ,M K (44) 

can be negative even if f(ps) is a regular positive function. For instance 
it is trivial to check that this happens for f(ps) = l/{p\ + m 2 ). There- 
fore De,hv{pe) is no more positive definite and clearly this is a potential 
source of problems for the perturbative expansion. 

2. The Wilsonian prescription of the double pole 

1 !™^T2 ( 45 ) 



\pl]w A-0p| + A 2 



is different from the CPV prescription (|42|). As a matter of fact the 
eigenvectors of the Euclidean propagator are explicitly non-negative 

°' p| + A 2 ' p| + A 2 ' pJ+A^ (46) 

and therefore one could argue that the Wilsonian prescription is better 
then the CPV. Unfortunately, the Wilsonian prescription has other 
problems that we will describe in the next section. 



4.2 The double pole problem in the Wilsonian formu- 
lation 

We found two major infrared problems in the Wilsonian formulation of the 
axial gauge: i) the Fourier transform of the propagator is divergent at A — > 0; 
ii) one loop Feynman diagrams which are infrared finite in covariant gauges 
becomes divergent in axial gauge for all configurations of momenta. In prin- 
ciple, this is not enough to prove the inconsistency of this gauge, since we 
should prove that at least one physical quantity is ill defined in the A — > 
limit. Indeed one could think that due to miracolous cancellations related 
to the Ward-identities (which are respected) physical quantities are indeed 
finite in the A — > limit as it happens for instance for the Wilson loop at 
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order 0(g 2 ) (see section 7). However this is not garanteed in general. More- 
over, even if this were true, the singularity of the A — > limit would forbid 
in practice any numerical application, since even extremely small numerical 
breaking of the Ward-identities would result in enormous (at the limit infi- 
nite) differences in the final result at A —>■ 0. This is a serious drawback of 
the axial gauge Wilsonian formulation. We stress that, even if not noticed in 
those references, the problems we discuss here also apply to the formulations 
in [16|, |17]] where they are probably even worse due to the wild breaking of 



gauge-symmetry. They were not recognized before simply because they do 
not affect the computation of the one-loop beta function, which is an univer- 
sal quantity not affected by the way the infrared cutoff in inserted and 
insensitive to the gauge fixing choice. 

The origin of all problems comes from the part of the propagator propor- 
tional to p^Pu, 

= PflPu n2 (47) 

A '"" (p 2 + A 2 )((p-n) 2 + n 2 A 2 ) 1 ; 

This quantity is a messy source of infrared divergencies due to the identity 
(in the sense of distributions) 

' ^° n 6(p-n) (48) 



(p ■ n) 2 + A 2 A 

This means that for any regular function f(p 3 ,A) such as /(0,0) ^ we 
have 

We will see in next subsections the disastrous consequences of this fact on 
the computation of various quantities. 

4.3 The x— space propagator. 

Here we show that even if the x— space propagator 



Da,^(x) = [ e- ip - x D A ^(p) 
J P 



(50) 
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is perfectly defined for any A ^ 0, the infrared limit A — > does not exist 
due to a strong infrared divergence coming from the double pole part of the 
propagator proportional to p^Pv (the double pole part proportional to A 2 n fJi n u 
and the single pole part are regular as A — > 0). This can be proved with a 
direct computation by using the result fl49f) . One obtains D^ i3 = D^ 3i = 
Av,33 = and 

A^O 1 [ ViVj -ip-x _ 1 lr r(3)/=\ A^O 



^i(-) = ^ /_ f^e"- = ^M (3) (*) A = U oo . (51) 



Therefore the Fourier transform of the propagator does not exist at A — > 0. 

The level of danger of this problem is unclear at this level, since it could 
not affect physical quantities. We will see in section 7 that the simplest 
physical quantity we can compute, i.e. the Wilson loop at order 0(g 2 ), only 
depends on the Fourier transform of the transverse part of the propagator, 

Dl,, v (x) = [ e~** f M (p) D x /(p) V(P), (52) 
Jp 

which is safe in the limit A — > 0. To show this point we have to prove that 
the D^ v> k{x) part of the propagator, 

Dk^{x) = DTixp^- (53) 

which is the only one which transverse part depends on the gauge vector 
is vanishing at A —>■ 0. A direct computation in Euclidean space gives 



DT(x E ) 



exp(ip E ■ x E ) 
„ L I'i: + A 2 



A 2 



pI + a 2 _ 



(54) 



By direct inspection using (48) one sees that this term is linearly vanishing 
with A; therefore the infrared limit of the transverse part of the propagator 
exists finite and is independent of the gauge vector n M : 

oU ^ s »( v _^)_l_ (55) 

The crucial point here is the fact that the singular term (^7|) does not 

contributes to the transverse part of the propagator. 
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Figure 1: Self-energy with two zero-momentum insertions. 



4.4 Computation of a one-loop integral 

The simplest perturbative quantity sensitive to the double pole problem is 
the one-loop quark self-energy. In this section we will see that the PpP v 
part of the propagator gives an infinite contribution to this quantity in the 
A — > limit. Unfortunately this is a gauge-dependent quantity therefore in 
principle we cannot positively conclude about the inconsistency of the axial 



gauge choice with the Wilsonian prescription ( |45| ) . However in practice this is 
a serious drawback since it forbids the direct computation of any perturbative 
quantity except the one-loop anomalous dimensions and the beta function 



In order to simplify the computation, by avoiding the inessential com- 
plications with the gamma matrices, we consider the self-energy of a scalar 
quark. Moreover, in order to avoid infrared divergences we take a quark 
mass m ^ and in order to avoid ultraviolet divergences we derive twice 
with respect to m 2 . In this way we obtain a quantity corresponding to the 
four-point vertex with two insertions at zero momentum shown in figure 1. 
This quantity, which we denote by 



f(p 2 ,pIA) 



'IE 



d 



dm 2 



{2p E + q E y{2pE + q E ) v 
[(q E +p E ) 2 + m 2 + A 2 } 



A) 



(56) 



is ultraviolet convergent and (in covariant gauges) infrared finite even at 
A = if p 2 E ^ — m 2 . Still, the A — > limit is singular in the massive axial 
gauge since the dominant contribution comes from the double pole part of 
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the propagator which, due to equation is infinity in the A — > limit: 



F(f,pl A) A =° --L j ^—±_ m{2 p + q y {2p + g y - oo. (57) 

Notice that this problem appears even in the Abelian theory. Moreover it 
is clear that the problem becomes worse in the non-Abelian case and at 
higher orders in perturbation theory. Still, in principle it is possible that 
cancellations allows to define unambigously physical quantities even in the 
A — > limit, since the various Feynman diagrams should add in such a way 
that the final result at A — > is finite and independent of the gauge fixing, 
i.e. the same that in covariant gauges. However, this is an highly non-trivial 
point and cannot be taken for granted at finite order in perturbation theory. 



5 The planar gauge 

For the reasons described above, we think that it is not convenient to in- 
sist on the axial gauge and from now on we switch to other gauge choices. 
The simplest examples where the difficulties of the axial gauge can be cir- 
cumvented (at least in the standard massless formalism and at order 0(g 4 ) 
in perturbation theory) are the planar gauge and the light-cone gauge. We 
consider first the planar gauge. 



5.1 Remarks on the planar gauge 

Properly speaking, even the standard massless planar gauge choice with the 
CPV prescription is not completely controled and the consistency of the 
perturbative expansion in this gauge is an open problem; nevertheless it will 
be studied in this section for sake of comparison with the axial gauge and 
the light-cone gauge choices. 

The reason of the difficulties can be traced back to the well known split- 



ting formula [13 1 



' ' + r / - ,m ) - n 2 5(nk)5(np) . (58) 



[nk] [n(p — k)\ np \ [nk] [n(p — k)\ 
Due to the presence of the delta function terms, in two-loops or higher orders 
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computations involving gluon vertices there appear ill-defined objects of kind 



1 kn s/tt 

\kn\ 5{kn > ~ (knf + e 2 (kn) 2 + e 2 (59) 

which have to be studied in detail. It is known that in the massless compu- 
tation of the 0(g 4 ) Wilson loop there is a precise way to manage these terms 



and the final result is consistent with the covariant gauge | I8[ . Neverthe- 
less, the situation at order 0(g 6 ) is unknown and in general the definition of 
two-loops Feynman diagrams is problematic. This is the reason why we will 
switch to the light-cone gauge in the next section. However it is interesting to 
see that even at the level of the planar gauge the divergences found before in 
one- loop computations are absent, essentially because the double pole term 



proportional to Pp!p v in (|47|) is absent. There is however a double pole in 



the D^ v part of the propagator, but since this term is multiplied by a A 2 
factor it is espected to vanish in the A — > limit. This is indeed the case for 
infrared safe quantities and it will be explicitly verified in the examples we 
consider here. We stress that this is by no means a trivial statement, since 
in the computation of particular infrared divergent quantities the term 
in the denominator could be not subleading with respect to the other terms 
and actually could also give the dominant contribution. We will see in sec- 
tion 6 that this is indeed the case in the light-cone gauge. Moreover we will 
see in section 7 that the D^ u term contributes to the Wilson loop and is 
suppressed only if the A — > limit is done properly, i.e. at finite T. 

5.2 The x— space propagator 

We prove here that the Fourier transform of the propagator is well defined and 
in the A — > limit reduces to the Fourier transform of the standard planar 
gauge propagator. This is obvious for the D^ v part of the propagator, 
therefore we have simply to prove that the Fourier transform of the D\^ u 
term is vanishing at A — > 0. A direct computation gives 

2 f exp(ip 3 x 3 + ip ■ x) A 2 ,. 2 _ 2 x n , Rn , 

lim A / = hm — ln(A x ) = 0. (60) 

A ->° Jp, P3 h] 2 G»i +P 2 + A 2 ) A-o 8tt 2 

This result can be obtained by first doing the p% integral using l/jp 2 ] = 
[p3]~ 1 d P3 ; the successive three-dimensional integral generates an infrared log- 
arithimic singularity at A — * which is killed by the A 2 prefactor. Therefore 
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we have solved the problem with the Fourier transform. However we stress 
again that in the computation of infrared sensitive quantities this term in 
general could gives an essential contribution to the final result. Furthermore, 
we will see in section 7 that this term is suppressed in the Wilson loop com- 
putation only if the A — ► and T — > oo limit are done in the correct order. 
Therefore the relation between the massive version of the planar gauge in 
the A — > limit and the standard planar gauge in general is delicate. 



5.3 Computation of a one-loop finite quantity 

Now we can study what happens in the computation of the quark self-energy 
with two zero-momentum insertions, i.e. the analogous of the quantity F 
defined in ([56]) . If we split the Euclidean planar gauge propagator in three 
terms (in this section for notational simplicity we write p for p#) 

D^{p, A) = — — - ISp, - 



p 2 + A 2 y [ p . n ] [p. n ] 2 j (61) 

= £>W(p, A ) + 2jW(p f A ) + Dgfa A ) 

we can split the computation in three different contributions: 

F(p 2 ,p 2 3 ,A) = F^\p 2 ,A) + F^(p 2 ,pI,A) + F^\p 2 ,pI,A) (62) 

with 

F^(n 2 \)= I ( 9 V ( 2 P + l) 2 ffi oN 

VP ' } J q \dm 2 ) {q 2 + A 2 ){{q + p) 2 + m 2 + A 2 )' 1 ' 

F V>(? p 2 A) = [ (^-) 2 -2(2p + q )-n(2p + q ). q 

J q \dm 2 ) [q-n}{q 2 + A 2 ){{q + p) 2 + m 2 + A 2 y 1 ' 

F (0 J A) = [ ( J_ V -A\2p.n + q .n) 2 

3 ' J q \dm 2 J [q-n] 2 (q 2 + A 2 ){(q + p) 2 + m 2 + A 2 )' 1 ' 

The first contribution is the same as in Feynman gauge and can be computed 
using the standard Feynman parametrization. It gives 

F (a )( , A) _ J_ f 1 dx (1 _ /(l+4x-, 2 ) + 2K(l-,)+A 2 ] 
{P ' ' 16tt 2 J 1 ] l(p 2 x + m 2 )(l-x)+A 2 } 2 

(66) 
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where the x— integral can be performed in terms of elementary functions 
(logarithms). In particular in the on-shell limit p 2 — > — m 2 , A — > the x ~ 1 
integration region dominates and we have 

F«V, A = 0) = + 0((i + m 2 f) (67) 

if the A — > limit is taken first. On the contrary, if the p 2 — > — m 2 limit is 
taken first, we have 

^V = _^ A) = _^_L + (£). (68) 

This is an explicit confirmation of the fact the the order of limits is a delicate 
question for infrared divergent quantities. 

The second contribution is typical of planar gauge. Its general expression 
is rather complicate, but there is a simplification if we restrict the analysis 
to the limits \p 2 \ ^> p 2 or \p 2 \ <C p 2 . The first limit \p 2 \ ^> p\ is trivial 
in the sense that the integral can be recast in a covariant-like form and its 
computation is no more difficult than the computation of the covariant term 
and gives a completely analogous contribution. In the second limit \p 2 \ <C p\ 
instead the peculiarities of the planar gauge are evident. In particular the 
quantity 

(pi A) = hm lF<V,p!) (69) 

can be computed with the method of double Feynman parametrization dis- 
cussed in appendix B. After some manipulation it reads 

with 

N = -2\j4(1 - x 2 ) 2 + p 2 3 [2(3x 2 - l)ql - q 2 (l - x 2 )/z] + qM + q 2 /z]] (71) 
and 

D = z z l 2 [q\ + q 2 + A 2 z + (pjx + m 2 )(l - x)} 3 . (72) 
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Now the integrations in q 3 ,q and z are straightforward (see appendix B) and 
all the complexity is confined into the x— integration. In the A — > limit 
there is a strong simplification and it can be exactly performed in terms of 
elementary functions. The final result is 

zr(b)f 2 a r\\ 1 p\/m 2 — 1 

V^i ln [( VpI + m 2 - y / p|)/( Vpf + y/pj + m 2 )} 
M 2 {p 2 + m 2 f/ 2 

The third contribution is characteristic of the Wilsonian formulation and 
must vanish in the A — > limit in order to recover the results of the standard 
massless planar gauge: 

Wm F {c \p E ) = 0. (74) 

This is immediate to see in the limit \p 2 \ <C p\. In the opposite limit \p 2 \ 3> p\ 
one could expect divergences at p 2 = — m 2 ; but it is immediate to see that 

F i - C \p 2 = -m 2 , A) ~ ^F^ a \p 2 = -m 2 , A) (75) 

therefore F^ c \p 2 = — m 2 ,A) is linearly vanishing at A — > 0. Thus the D\^ v 
term can always be neglected, at least in this kind of one-loop computations. 
We will see that this is not always the case for the light-cone gauge. 



6 Safeness of the light-cone gauge 

For the point of view of the consistency of the perturbative expansion the 
light-cone gauge is expected to be the safest choice. There are various reasons 
for this expectation, which we will discuss now. 

The first reason is the fact that, contrary to the planar gauge, the Eu- 
clidean propagator is (semi) positive definite. This can be immediately proved 
by solving the eigenvalue equation 

det (\5Z-DZ(p E ,n E ; A)) = (76) 

which after explicit computation reads 

i ^ 2 



~ p|h-a^J =° (77) 
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with solutions 



Ai — A2 — —5 — ; — T-r, A3 — A4 — 0. (78) 

p E + A 2 

From this computation one sees that the eigenvalues Aj = \ijpEi t^e) are non- 
negative. In addition, they depend only on the Lorentz-invariant combination 
p E and not on tie- We emphasize that this is not true in non-covariant 
gauges others than the light-cone one. Therefore one expects a somehow less 
prononced Lorentz breaking in the light-cone gauge. 

The second reason is that within the ML prescription the splitting formula 



( |58| ) holds without the delta function term: as a consequence there are no ill 
defined terms in higher orders of perturbation theory as it happens with the 
CPV prescription in axial and planar gauges. 

The third reason is that the massless limit of the light-cone gauge is ex- 
pected to be less singular than in others gauges. This expectation comes from 
the that the Fourier transform of the D\^ v term is quadratically vanishing 
as A 2 — > 0: this should be contrasted with the axial gauge case where the 
suppression of the D\ )fll/ term in only linear in A (see the discussion after Eq. 
(|54|)) and the planar gauge where this term vanishes as A 2 times a logarithm. 
We will discuss in detail this latter point just below. 

6.1 The x— space propagator 

In order to explicitly compute the Fourier transform of the D^^ u {x) term we 
have to study the Euclidean integral 

DT(x E ) = -A 2 / (79) 

Jg E {qE-n E ) 2 {q E + A 2 ) 

We proceed as follow. First, we introduce the angles 9 and 9' such as 

<?i = cos 9', q2 = q± sin 9', q% = q\\ sin 9, q^ = q\\ cos 9 (80) 

and we take x = (x±, 0, x\\, 0) (this is not restrictive); then the Fourier trans- 
form can be written 

/n Ar, API 1 An,, AO o~ 2i9— iq\\x\\ sinfi— iq±xj_ sin#' 
q ± dq ± d9 dq\\dV e n II 
^TJ2^ ,,,^ + gi + A 2 ) • (81) 
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D™(xe) can be computed by first performing the angular integrations and 
then the momentum integrations in q\\ and q±. The angular integrations can 
be done by using the following representation of Bessel functions 



JJz) 



i r 2n 

— / de e - ine - lzsint 



for n = and n = 2 respectively. In this way we obtain 

A 2 r , f°° dq\\ J 2 (q\\x ll )Jo(q±x ± ) 



DT(x E ) 



(27T) 



o q\\ q\ + q\ + A 2 



(82) 



(83) 



Since at small z the Bessel function J-iiz) — g- 2 is quadratically vanishing 
whereas at large z is exponentially suppressed, we see that the gy— integral 
is infrared and ultraviolet finite. It can be exactly computed with the result 



DT(x E ) 



A 2 



u 2 , - 4a% 2 + 2a; 2 K 2 (uj ± x\\) 
dq±q±— S ^Jo(g±x ± ), (84) 



_L 



with a;^ = g 2 + A 2 . Using the series expansion of the Bessel functions 



Ko(z) 



2z 2 



+ 0(z 2 log z 2 ), J (z) 
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one sees that the q±— integral is also finite, even in the A — > limit where it 
can be computed exactly and it gives 



DT{x E ) 



A 2 



16tt 5 



4io g 4-i + o(AVLA 2 zjT 

x \\ 



xi 



(85) 



where x 2 E = x 2 ^ + x 2 ± . Thus, Da^ u (x) is quadratically vanishing in the A — > 
limit. We also see that there could be problems in the limit x\\/x± — > oo which 
is relevant for the Wilson loop computation. This point will be discussed in 
detail in section 7. 



6.2 The structure of (soft) infrared divergences 

It is possible to repeat the computation of section [5]3] in the light-cone gauge. 
The covariant contribution (p, A) is obviously the same, whereas the 
contributions F^ b \p,A) and F^ c \p, A) are different from the analogous one 
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in planar gauge. Still, one can use the double Feynman parametrization 
method described in appendix B. In the limit p\/p 2 — ► the contribution 
can be rewritten in the form 

F»(ri a) = f t dx f dz (J-) 2 jMl^flj (86) 



3||.9± •'O ■/1-a; 



with 

iV(gf,?i,p|,x,z) =-8gfpf + AiV( ? J,gi,pf,x,z) (87) 

and 

D = + q\ + A 2 Z + + m 2 )(l - x)] 3 . (88) 

Actually, we exactly computed the expression AN(q^,q\,p^,x, z) in the nu- 
merator and it is possible to perform exactly the integrals in the general case, 
but the resulting expressions are lenghty and they will not be reported here. 
In order to analyze the on-shell divergences at p 2 ~ — m 2 the only relevant 
term is the fist one displayed in the right hand side of fl87D . The integrations 
in q\\ and q± are straighforward and they give, after neglecting the AN term, 

F^(p 2 A)^--L f dx f -j- /f 1 "*?' ( 89 ) 

v 11 ' Ait 2 J J x _ x z [A 2 z+ (p 2 x + m 2 )(l -x)] 2 

If we take the A — > limit first we obtain 

nW , 2 . nxP 2 — m 2 1 hi[(pj + m 2 )/p 2 ] / i \ 

and this is dominant with respect to the divergence l/(p 2 + m 2 ) coming from 
the covariant contribution F^ a \ If we take the p 2 — > — m 2 limit first we 
obtain 

^=^=ps +o O> (9i) 

which is still dominant with respect to the linear divergence coming from the 
covariant contribution F^ a \ 
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Now we can study in an analogous way the contribution of the D\^ v 
part of the propagator, i.e. the term. It is obvious that this term is 
quadratically vanishing at A — ► if p 2 7^ — m 2 . Nevertheless, if p 2 = —m 2 , 
this term could be divergent as A — > 0. This is the case indeed. Actually 
it is easy to understand that there is no divergence in the limit [pjj| -C p 2 ^, 
using the same argument as in equation (|75|) . However in the opposite limit 
\P\\\ ^ P 2 ± ^ ne infrared divergence is enhanced and actually dominates with 
respect to the covariant contribution. To see this point we write down the 
explicit expression of in the near on-shell region, obtained with the 
method of double Feynman parametrization. It reads 

F^ipl A) =f fdx I" dz flVj*M!il_ (92) 
"' J^Jo Ji-x \dm 2 ) D( q 2 ,ql,p 2 ,m\A^x,z) 

with 

N(ql q 2 x ,pl x, z) = 24(1 - x) 2 p\A 2 + AiV(gg, p% x, z) (93) 

and 

D = z\q\ + ql + A 2 z + (p\x + m 2 )(l - x)] 4 . (94) 

In the term AN(q 2 , q\,p 2 , x, z) in the numerator we collected all the contri- 
butions which are sub leading in the on-shell limit. After double derivation 
with respect to the mass m 2 and momentum integration we obtain 

1 fi dz A 2 K(1 - ^) 4 



x_ x z [A 2 z+ (pjx + m 2 )(l -x)} 4 ' 

(95) 

In order to study the A — > limit at p 2 = —m 2 it is convenient do rescale 
the z variable as 

f 1 dz f 1/{1 - x) dz , . 

z = z(l-x), / — = / — (96) 



l-x 



Z 



We see that for x — > 1 the upper limit of the z— integral tends to infinity, 
therefore the z— integration simplifies. Then the x— integration is done taking 
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in account that the x ~ 1 integration region dominates in the near on-shell 
region; finally one obtains 



F^ c \p\ = -m 2 ,A) 



1 m 2 



(97) 



6tt 2 A 4 ' 



It is important to notice that even if the Da term dominates the soft di- 
vergence at m 7^ 0, it can be always neglected in the analysis of hard or 
collinear divergences, i.e. if m — 0. In this case in fact the contribution 
from the -Da term is quadratically vanishing. This point can be understood 
with a dimensional argument and explicitly checked by using the double 
Feynman parametrization formula and by noticing that at m 2 = it is the 
x ~ region of integration which dominates, not the x ~ 1 region. Then the 
z— integration can be done by using 



and the x— integration by using the tricks reported in appendix B. With 
similar techniques we checked that the gluon self-energy integral in the A — > 
limit goes smoothly to the standard massless integral for off-shell Euclidean 



The Wilson loop is the simplest physical quantity where the effects and the 
problems of the infrared regularization can be studied. In particular our 
scope here is to study the subtilities of the A — > limit and to test how 
the essential property of the gauge-invariance of the Wilson loop is recovered 
when the infrared cutoff is removed. In concrete in this section we compute 
the Wilson up to order 0(g 2 ) in perturbation theory. This is enough for 
elucidating various important features of massive axial, planar and light- 
cone gauges and it is a first step versus a more comprehensive computation 
at order 0(g 4 ) in perturbation theory. 

For definiteness, we shall consider a rectangular Wilson loop T LT of size 
2L x 2T, with T ^> L. We shall work in Euclidean space with coordinates 
xi,X2,X3, X4 and we shall take the loop in the plane X2X3, as shown in figure 
2. We should notice that T denotes a lenght in the spatial direction x 3 , i.e. 




(98) 



momenta p 2 7^ — p\. 



7 The Wilson loop test 
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Figure 2: Rectangular Wilson loop of size 2L x 2T. 

the direction of the gauge vector = (0,0, 1,0), and not in the temporal 
direction X4 = ixo. Therefore apparently this loop in not related to the 
interquark potential. Nevertheless if the theory is consistent all directions 
in the Euclidean space must be physically equivalent at A — > 0, therefore 
the final result must be the same that for the loop corresponding to the 
interquark potential where the side T is in the temporal direction. 
Formally the Wilson loop is defined by 

W rLT = ^< TrPexp (ig jf A^daf) > (99) 

where P denotes the path ordering on the loop Tlt, Tr is the trace in the 
fundamental representation of SU(N C ) and the average is evaluated via a 
perturbative expansion of the Euclidean functional integral. In our case Y LT 
can be split in four pieces 

r LT = r« + r( 2 ) + r( 3 ) + r( 4 \ (100) 
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parametrized as 

Wis) 



sT 



r (3) ( s ) = sT 



( o \ 



+1 

V o / 
/ o \ 



-1 

V o / 



r(2) ( . 



r (4) ( 



sL 



sL 



( \ 

-1 


V o / 
/ o \ 

+i 



V o / 



(101) 



where the parameter s lives in the interval [—1,1]. Equivalently, we can 
parametrize 



r«(*) = *U if (i) 
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Tf>{t) = -t6p, r«(i) = *U (102) 



with t G [— T, T] and I G [—L,L]. The Wilson loop can be easily computed 
by expanding in powers of the coupling constant, 



LT t LT 



9 -W^ + i-W® + Otf) 
2 3! 1 lt vi( ' 



(103) 



By using the properties of the generators in the fundamental representation 
of SU(N C ), Tr 1 = N c , Tr T a = and Tr(T a T b ) = |<5 afe , from (|99|) one obtains 
the explicit expression 



AN, 



dx^dy" < A^x) ■ A u (y) >W +OQ7 4 ). 



(104) 



c JT 



Notice that higher orders corrections begins at order 0(g 4 ), not 0(g 3 ). The 
expectation value < A^(x) ■ A v (y) > ( -°\ computed at zero order in the cou- 
pling constant, coincides with the Fourier transform of the free propagator, 
therefore a a more explicit expression for the 0(g 2 ) contribution is 



wj? ] = -C F 

1 LT F 



1 dr^ &v v 

ds 2 - —D^(T(s 1 )-T{s 2 )), 



i 



dsi ds2 



(105) 



where the relation 5 aa = 2N C C F with C F = (iV c 2 - l)/(2iV c ) has been used. 
Specializing to the rectangular Wilson loop in figure 2 we have in general six- 
teen contribution to (|105|) ; however many terms gives the same contribution 
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i LT 

2C F 



due to the symmetries D^x) = D VfM (x) and D^ v (x) = D^ u (—x) and finally 
one is left with the explicit expression 

/ dti / dt 2 [1)33(0, 2L,*x + t 2 ,0) + ^33(0, 0,^-^,0)] 

J-T </-T 

dZi y^d/ 2 [£> 22 (0, Zi + Z 2 ,2T,0) + ^22(0, Zi-Z 2 , 0,0)] 

+4 y dt y dl D 23 (0,L + l,T -t,0). 

(106) 

This formula can be further simplified in the T — > 00 limit since various 
contributions are subleading. Moreover, due to the identity (which holds 
since T LT is a closed path) 

/ dafdtfD^x - y) = f dx»dy v Dl v {x - y), (107) 

J Tlx J^lt 

actually only the transverse part of the propagator contributes. At zero mass 
Dj^ u (p) is the same in all gauges and this is the reason why the final result 
is gauge- independent; however at A 7^ there is a dependence on the gauge 
vector n M coming from the D Alxv part of the propagator. 

A naive expectation would suggest that D\,ni> being proportional to A 2 
gives a vanishing contribution to the Wilson loop. This is in general should 
not be taken for granted. Nevertheless, suppose for a moment that this naive 
espectation is correct (this point will be analyzed in detail in next section). 
Then one can consider only the D^ v part of the propagator or even only the 
transverse part D^ u . This object is independent of the gauge- fixing vector 
n M and therefore the Wilson loop is the same as in covariant gauges. In 
particular, we can effectively replace D^ v with 5 fJil/ /(p 2 + A 2 ) = S^D A (p) 
because they have the same transverse part. If we divide by 2T in order to 
eliminate the contributions subleading at T — > 00 we obtain for the 0(g 2 ) 
analogous of the interquark potential 

(2) 9 2 N c Wj 2) 
vi 2) (2L) = -lim y C 4t Tlt (108) 

the expression 

V1 2) (2L) = 9 —^ J dt x J ^dt 2 D A (t 1 + t 2 ,2L)+c (109) 
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where c is the infinite renormalization constant 



T^oo 2T 



T pT 



a N C 

lim ^ F I dtx I dt 2 D A (t 1 -t 2 ,0). (110) 

-T J-T 



This is the well known ultraviolet divergence of point-like charges and can 
be eliminated by fixing the potential to be zero at L — > oo. With simple 
manipulations involving the representation of the delta function 

2n5(p 3 ) = lim / dx 3 e ipsX3 = lim 2sm ^3 T ) (m) 

T^oo J_ T T^oo p 3 



and the formula 

(2sinp 3 T) 2 = 2T 

T^oo pj T^oo KF6h ^ 1 

from (|109|) one obtains 

if > (2i) = c /n cCf I = fNJ3 ~*g£±. (113) 

This is the formula of the screened Coulomb potential of two colored charges 
at distance 2L and reduces to the standard one when the infrared cutoff is 
removed. We still stress that this formula is correct as far as we can neglect 
the contribution from the D\ term in the limit A — > 0. 



7.1 Vanishing of the term 

Now we prove that the D\ term of the propagator gives a vanishing contri- 
bution to the Wilson loop. This is a non-trivial point since in this quantity 
there is a problem of commutativity between the limits A — > and T —>■ oo. 
This can be seen from the explicit expression of the D\ contribution to 
the interquark potential, which up to an unessential infinite constant reads, 
working for definiteness in the planar gauge, 



V1 2) (2L) = J* dt x J* dt 2 D 33 (0, 2L, *! + t 2 , 0) 

g 2 N c C F f e ip22L A 2 (2sinp 3 T) 2 



(114) 



2T J PE [p 3 ] 2 (p 2 3 +P 2 + A 2 ) p\ 
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We see that we cannot blindly use the formula ( |112| ) and compute the integral 
directly at T — > oo. We are instead forced to consider finite T and to make 
some subtle observation. The idea is that the p%— integral is dominated by 
the P3 ~ region and that actually it reduces to an integral of the kind 



UT) = L^r^ (115) 

with a = 1. Strictly speaking this is a divergent integral; nevertheless it can 
be defined in the region —1/2 < a < 1/2 where its value is 

I a (T) = -r(-l - 2a) sin(a7r)(2T) 1+2Q . (116) 

7T 

After analytic continuation to a = 1 we obtain 

J 1 (T) = ~T 3 . (117) 



A simpler way to derive this formula is to derive equation (|117|) three times 
with respect to T: then the sinus representation of the delta function is 
recovered and the integral is easily computed. In other words we can use the 
formula, which has to be interpreted in the sense of distributions, 

(2sinp 3 T) 2 4 » r/ . , . 

lim 1 \\ } = lim --T 3 • 2n5(p 3 ). 118 

Then the three-dimensional integration is as in ( |113j) and the final result is 

^ = 2g*N c C F exp(-A.2L) 
A v ; 3 4tt-2L v ; 

We see that this contribution is quadratically vanishing in the limit A — » 0, 
at finite T. The T — > oo limit cannot be taken before the A — > limit. This 
is the crucial point of our analysis. 

For what concerns the situation with other gauges, we observe that in the 



light-cone gauge we have to compute exactly the same integral as in (\1 14ft 
and therefore we have the same result. On the contrary, in the massive axial 
gauge we have to replace a factor l/[p3] 2 with a factor (pjj + A 2 ) -1 and using 
(|4lf) we see that the final result is linearly vanishing with AT. This means 
that the 0(g 2 ) Wilson loop test works for any gauge choice: but this is hardly 
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a surprise since the 0(g 2 ) computation corresponds to the computation in 
an Abelian theory. It is well known that problems begin in the non-Abelian 
case and at higher orders in perturbation theory, starting from order 0(g A ), 
therefore we cannot conclude nothing about the final consistency of these 
Wilsonian gauge choices at this level. Nevertheless we think that this anal- 
ysis of the free case is very instructive and allows to learn a lot about the 
possible origin of problems. In particular we have learned that in an 0(g A ) 
computation the only part of the propagator we have to control is the D\ 
part, which is quite simple and surely can be studied with a relatively little 
analitical effort. In other Wilsonian formulations of non-covariant gauges 



based on infrared cutoffs more complicate than a mass-like term [|16| , |17 
the analysis is technically much more cumbersome but physically equivalent 
in what concerns the final results in the A — > limit. The drawback of 
the generic cutoff is that the fine-tuning problem must be solved in order 
to fix the correct boundary conditions (renormalization prescriptions) such 
to have an infrared limit consistent with the gauge-symmetry; this difficult 
problem is avoided with the mass cutoff since with this choice the theory 
is Ward-identities-consistent to all scales. Finally we would stress the fact 
the planar gauge and in particular the light-cone gauge are expected to be 
much more regular than the axial gauge in the zero mass limit: this expec- 
tation is reflected in the present computation by the fact that the D\ term 
is quadratically suppressed in both planar and light-cone gauge and merely 
linearly suppressed in axial gauge. 



8 Conclusions 

In we pointed out that in algebraic non-covariant gauges it is possible to 
build up a Wilsonian formulation of gauge theories consistent with the Ward 
identities provided that the infrared cutoff is introduced as a formal "mass" 
term. However, we stressed that this infrared cutoff cannot be physically 
interpreted and must be removed in order to recover the essential property 
of gauge-independence of physical quantities. In this paper we have inves- 
tigated the properties of the singular limit A — ► by explicitly computing 
various quantities. We have seen in general that the pure axial gauge choice 
is problematic since the Fourier transform of the propagator and even the 
simplest Feynman diagrams which are finite in covariant gauges are instead 
divergent at A — > 0, for any configuration of momenta. Moreover, we argued 
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that other gauge choices, namely the planar gauge and the light-cone gauge 
are in a much better shape. We have however seen that the structure of 
infrared divergences is quite subtle even in these cases and that there are 
gauge-dependent quantities such as for on-shell configurations of momenta 
the contribution from the gauge- dependent term of the propagator D\^ v is 
not only non-vanishing, but even it is dominant with respect to the standard 
contribution. Therefore the A — > limit for these quantities is delicate. How- 
ever these quantities are unphysical and this fact should not be considered 
as suggesting an inconsistency of the theory. The only way to prove the con- 
sistency or the inconsistency of the approach is by considering true physical 
quantities like the interquark potential as obtained from a Wilson loop of 
size 2L x 2T with T — > oo. We have seen that the limits A — > and T — > oo 
must be studied with great care, since they do not commute. In particular 
the crucial property of the exponentiation of the Wilson loop is recovered 
only if the A — > limit is taken before the T — > oo limit. We have seen that 
the planar gauge and the light-cone gauge appear to be much more regular 
than the axial gauge in the infrared limit. This is explicitly realized in the 
Wilson loop 0(g 2 ) computation by the fact that the gauge dependent term 
is quadratically suppressed at small AT whereas in axial gauge this term is 
only linearly suppressed. We notice that in covariant gauges, by accident, 
the order of limits is not crucial at order 0(g 2 ), then this feature has not 
been recognized previously; nevertheless it is manifest in an 0(g 4 ) compu- 
tation. The reason is that the infrared cutoff breaks the BRST- invar iance 
and therefore the cancellation of the so called non-Abelian contributions to 
the Wilson loop no more works at A ^ 0. Therefore the interquark potential 
cannot be perturbatively defined at A ^ and T — > oo, but only at finite 
T . Finally we notice that the physical interquark potential, defined exactly 
at A = 0, is independent of the cutoff function choice when it is computed 
order by order in perturbation theory. 

For the future, we plan to study in detail what happens in the 0(g 4 ) 
computation. Actually, at finite T, we expect that in planar and light-cone 
gauges the computation reduces smoothly to the standard massless compu- 
tation, which is consistent, whereas the axial gauge should be very delicate 
and possibly inconsistent. Still, we should notice that even in the cases where 
the 0(g 4 ) Wilson loop test fails, in principle this only indicate a failure of the 
perturbative expansion and not necessarily of the full theory. Nevertheless 
one would feel much more confortable with an approach admitting a pertur- 
bative expansion. In such a perspective, the most promising possibility seems 
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to be the light-cone gauge which certainly deserves additional investigations. 



Note Added. During the completion of this work we received a com- 
munication from R. Soldati and A. Panza who explicitly computed the 
Wilson loop at order 0(g 4 ) in the massive axial gauge case and proved that 
the A — > limit is singular. Therefore the axial gauge choice seems to be 
definitely pathological even in the Wilson renormalization group approach, 
at least at the perturbative level. 
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A Gauge-dependence of the pressure 

In this appendix we provide a very simple and illuminating example which 
illustrates the gauge-dependence problem in presence of a non-zero mass 
cutoff A. 

Consider a free gas of photons in thermal equilibrium at temperature 
T = 1/(3 in a box of volume V. We can compute, as a typical quantity 
directly related to the partition function, the pressure of this gas in presence 
of the mass cutoff A. We obtain different result in different gauges, but 
all these results collapse to the correct physical result in the physical limit 
A — > 0. This is trivial example, since the theory is free, nevertheless we 
believe it is very instructive. We recall that the pressure is a typical quantity 
which cannot be computed in perturbation theory in thermal field theory 
and where the non-perturbative powerfulness of the Wilson renormalization 
group approach could give an alternative way of facing the problem. In this 
sense it is an interesting quantity. 

We begin the computation by recalling some elementary facts (see for 
instance [0]). In quantum field theory the thermodynamic pressure 



p(P,V) = ^lnZp V (120) 
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is defined in terms of the partition function 



Zpv = J [#] exp (- J dr J d 3 x C E (<f>, d^- A )) (121) 

where Ce{4>, d^; Ao) is the bare Euclidean lagrangian of the theory, which 
is a function of the bare parameters depending on the ultraviolet cutoff A . 
For instance for a free scalar field of mass A we have, 

C E {<t>, d^- Ao) = Q^w) + l -A 2 <\? + c 4 (A, Ao), (122) 

where the notation Aq reminds that there is an ultraviolet cutoff 

inserted in the propagator in momentum space and the term c 4 (A,A ) = 
Aq[c 4 + C4A 2 /Aq + 0(A 4 /Aq)] is a vacuum energy counterterm of dimension 
four, which will be fixed later on by imposing the normalization condition 

lim In Zp V = 0, (123) 

i.e. the partition function is fixed to be 1 at zero temperature. It is interesting 
to notice that in order to impose this normalization the ultraviolet regular- 
ization is needed even if the theory is free. Using the gaussian integration 
formula we obtain 

Z pv = Afdet(-d 2 E + A 2 )' 1 ' 2 exp(-c 4 (A, A )(3V) (124) 

where J\f is a temperature independent normalization factor to be fixed later 
whereas the determinant of the operator — d 2 E + A 2 (acting on functions pe- 
riodic in the imaginary time) is defined as 



det(-9| + A 2 )^ 1 / 2 =exp 

pv 



-ilVln(-di + A 2 ) 



(125) 



with 



oo „ 

Tr?vH-d 2 E + A 2 ) = V J2 /m[(?4) 2 +p 2 + A 2 ], pi = 2vmT. (126) 

n=-oo J P 
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Using the stardard summation formula^ 

oo 

H(Pn) 2 + X 2 ] = -2 log[l + n{X)\ + (3X, (127) 

n=— oo 

where n(X) = (exp(/3X) — is the Bose-Einstein distribution function, 
one obtains 

In Zp v = V { m[l + n(w p (A))] - \(3V [ ° oo p (A) - (3Vc 4 + In A/", (128) 
■Jp ^ Jp 

with o; p (A) = a/p 2 + A 2 . The first term vanishes at /3 — ► oo. Imposing 
the normalization prescription ( |123|) fixes the counterterm c^A, Aq) and the 
normalization factor M to be 

c 4 (A,A ) = ~j °^(A), JV= 1, (129) 

therefore the pressure reduces to the well known expression 

P(T, A) =T | ln[l + n{u p )\ A =° ^T 4 . (130) 

Now we can repeat the same computation in the gauge theory case which is 
less trivial. We will consider the general linear gauge (|T6|) with £2 7^ 0. Using 
the Gaussian integration formula we obtain 

ZpviL, 6) = ATdet £ 2 - 1/2 det d^L„ det D X J U 2 (131) 

p V pv 

where det^ 1 ^ 2 comes from the integration of auxiliary fields, detpy d^L^ 

1 12 

from the integration of ghost fields and det^-D^i, from the integration of 
gauge fields. With a lengthy but straighforward computation, one obtains 
the determinant of the Euclidean propagator in the general class of linear 
gauges as 

det D A ^( P ) = (p 2 + A 2)2 (p . L) 2 + A 2[/ 2 2 + 6(p2 + A 2)] • ( 132 ) 



7 In order to use equation ( |127|) we have to remove the ultraviolet cutoff in the energy 
variable 734, otherwise the summation reduces to a finite sum with TV ~ Ao/T terms. In 
general the summation formula is defined up to a possibly divergent constant which can 
be reabsorved in the vacuum energy counterterm C4(A, Aq). 
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Suppose for a moment that A = 0: in this case we see that the partition 
function is gauge-independent since the £ 2 term in the photon propagator de- 
terminant is cancelled by the contribution coming from the integration of the 
auxiliary fields whereas the L— dependent part is cancelled by an analogous 
contribution from the ghost propagation, which in turn is fixed by requiring 
BRST-invariance of the total action. However, as can be explicitly seen, this 
cancellation mechanism of the gauge-dependence does not work in presence 
of an infrared regulator which breaks BRST-invariance. Therefore at finite 
A one obtains an unphysical dependence on the gauge-fixing parameters. In 
particular in the limit £2 - > i.e. strictly imposing the condition L^A^ = 
one obtains for the logarithm of the partition function the explicit expression 

In Z pv = 2Trp V ln(-d% + A 2 )~ 1/2 + Tr^ ln((<9 • L) 2 + L 2 A 2 )~ 1/2 
- TV ln((<9 • L) 2 y 1/2 ~ PVc A (A, A ) + \nAf . 

The first term is the expected one, corresponding to the pressure of two 
bosonic degrees of freedom; however there are also terms which are explicitly 
L— dependent at A 7^ and therefore unphysical; the L— dependence only 
cancels at A — > where one recovers the correct result 

p(T;L) A =°^T 4 VI. (134) 

This is the generic situation, however it is interesting to study what hap- 
pens in specific gauge choices since the gauge- dependence of the pressure 
can be less prononced than expected. For instance in the light-cone gauge 
L 2 = n 2 = we see that the n^— dependent terms cancel even at A 7^ 0. Ac- 
tually the n M — dependence cancels in the large class of non-covariant gauges 
such as L M (p) does not depend on the p 4 variable (for instance this is the case 
for the planar gauge and the Coulomb gauge): in this case the L— dependent 
terms in In Zpy have the form j3Vf(L, A, A ) and therefore can be reabsorved 
in the in the vacuum energy counterterm c^A, A ). In other words, they are 
eliminated by the normalization condition ( |123| ) and do not contribute to the 
pressure even at A 7^ 0. However the dependence on the quantization direc- 
tion should be expected in a two-loop computation; moreover it is evident in 
other observables, for instance in the Wilson loop computation, therefore in 
any case the infrared cutoff A cannot be physically interpreted as a true mass 
term, consistently with the standard lore that the only way to give a physical 
mass to a non-Abelian theory is via the Higgs mechanism. Nevertheless from 
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this simple example one could extrapolate the conjecture (to be checked case 
by case perturbatively with higher order computations or non-perturbatively 
with a numerical analysis) that for particular observables there are classes 
of gauges more regular that others, in which the gauge- dependence is very 
mild even for non-zero A. Clearly, this is a very interesting point in the 
spirit of phenomenological numerical analysis and should be investigated in 
the future. 



B One-loop integrals 

In this appendix we give some generalities on the computation of Euclidean 
one-loop integrals in planar and light-cone gauges. We begin by fixing our 
notations on integrals: for axial and planar gauges we define 

f=f£h, / = /=/t^. (135) 



„, ., W J„ J 2n' J, J W 
with q = (q°, q 1 , q 2 ), whereas for the light-cone gauge we define 

d?q± f f d 2 q\\ 



(27T)" J n J (2tt)" 



(136) 



with q± = (g 1 ,? 2 ), q\\ = (q 3 ,q )- Euclidean vectors are obtained after Wick 
rotation p 4 = ip ; we shall use the notations 

Pe = (p, ipo), Qe = (q, iqo), PeQe = V V E Ve = ~ 9 = ~PQ ( 137 ) 
and 

P\\ = Pi, P± = P 2 (planar gauge) (138) 



P\\ = P\ + Pi p\=p\+ p\ (light-cone gauge) (139) 

In this appendix we will always work in Euclidean space even if for sake of 
notational convenience the index E will be neglected. 

For semplicity we will restrict our remarks to the computation of the 
one-loop self-energy of a scalar quark. In general this is a rather complicate 
function F = F(p\\,p±,A), but we can give analytical estimations in the two 
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limits |p|| <C p\ and \pV\ ^> p\. In the first case by putting p\\ = we 
see that the spurious terms • n] cancel and then the one- loop integral 
can be put in an explicitly covariant form: therefore it can be computed 
through the usual Feynman parametrization. The case [pjj| 3> p\ instead 
is more cumbersome, nevertheless it can be easily implemented with a sym- 
bolic manipolation package. As a matter of fact, we prepared a set of routines 
based on the double Feynman parametrization method allowing to compute 
analytically all the integrals which are encountered in one-loop self-energy 
diagrams, including finite parts, and we checked that the A — > limit re- 
produces the known results of the standard approach. However, for sake of 
brevity, here we simply sketch the analysis for the simplest examples: the 
generalization to more complicate cases is straighforward. 

The general form of the Feynman integral one encounters in the evaluation 
of the quark self-energy (or its derivatives with respect to the mass) is given 
by the expression 

F((p ■ n? p 2 A) = [ N(q,p,n) 

The integral ( |140| ) will be explicitly evaluated both in planar gauge and light- 
cone gauge in next subsections. 



B.l Planar gauge integrals 

Consider first the planar gauge case where q ■ n = q%. The most conve- 
nient technical tool we found to manage this kind of integrals is the double 
Feynman parametrization which consists in using the identities 

' f d v r , . ,S . a (141) 



q 3 q 2 + A 2 J [q 2 + y(q 2 + A 2 )] 
and 

N(q,p,n) = 

(q 2 + A 2 ) a {{q + p) 2 + m 2 + A 2 ) b 

1 N(q — p(l — x),p, n)x a_1 (l — x) b ~ l 



B{a,b)J q J [q 2 + Ax + (q 2 +p 2 x + m 2 + A 2 )(l-x)] a+b . 



(142) 



Then it is convenient to introduce the variable z = (y — l)x + 1 lying in the 
interval 1 — x < z < 1 and to rescale q = z x ^ 2 q\ in this way equation ( 140 ) 
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can be rewritten in the form 



LJo h-* z^ 2 [q 2 + q 2 + (p 2 x + m 2 )(l-x)+A 2 z}^ ^ 

where we have introduced the translated and rescaled numerator 

[q 3 N(q,p, n)] tr = [q 3 N(q,p, n^J-Mi-,) ■ (144) 

After symmetrization in q 3 and three-dimensional angular average in the 
numerator, the momentum integrals can be performed by using the general 
formula 

(ql) Ml (q 2 ) M2 _ Bi(Mi,M 2 ,N) 



> , r ,. q (?f + <f + A ) N 8ir 3 A N - M i- M *- 2 

where Bi(Mi, M 2 , N) denotes the product of beta functions 

B x = B(N -M 1 -M 2 -2,M 2 + 3/2)B(N - M x - 1/2, M x + 1/2). (146) 

The integral in z is trivial at A = and a little complicate at A ^ but 
still expressible in terms of elementary functions; the integral in x instead 
in non-trivial. Nevertheless one can explicitly check that it is finite and in 
general can be expressed in terms of special functions. In the A — > limit 
the analysis strongly simplifies and one obtain the esplicit result reported in 
the text. 

With a simple generalization of this method one can compute Feynman 
integrals where the double pole l/[q ■ n} 2 appears, without encountering any 
problem. Moreover, the extension of this method to the integrals appearing 
in the gluon self-energy computation is straightforward. 

B.2 Light-cone gauge integrals 

Consider now the light-cone gauge case where q-n = iq^ + qz- Still the double 
Feynman parametrization can be used, it is enough to replace identity (|141|) 
with 

d V TZ2 I 12 I ,./„2 i A2M9 - ( 147 ) 



m + qsq 2 + A 2 J [qj + q 2 3 + y{q 2 ± + A 2 )}' 
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Using ( |142| ) , introducing the variable z — (y — l)x + 1 and rescaling q± = 
z 1 / 2 q± the Feynman integral (|140|) can be rewritten in the form 

'\ T f 1 dz [Hg4 + g3)^(g,P,n)Ul-x)° 



<?_L •'O </ 1-x 



z[q 2 + q\ + (p 2 x + m 2 )(l — x) + A 2 ^| 



where we have defined the translated and rescaled numerator as 

[(-% + g 3 )A%p,n)] tr = [(-ig 4 + qz)N{q,p,n)]f^ Z n{1 _ x y (149) 
Now we can perform the angular average in two dimensions 

1 r2n 



< f(q\\ -P\\) >n 2 = I d9 f(q\\P\\ cos 6), (150) 



o 



by using the general formulae 



< (si • Pi) 2 " >n 2 = J^r^ 2 ]) € pf> < (ffn • Pl) 2n+1 >n 2 = 0. (151) 

In particular 

1 3 

< -P\\) 2 >n 2 = 2 gy f>jj> < • P||) 4 >n 2 = g Q\\ P\\- (152) 



The integrals in gy, can be performed by using the general formula 
/■ (9 2 ) Ml (?l) M2 5 2 (M 1 ,M 2 ,iV) 



(g 2 + g 2 + A)* i 67r 2 A JV-M 1 -M 2 -2 ( 153 ) 



where 

B 2 = B(N -M 1 -M 2 -2,M 2 + l)B(N - M 1 - 1, M x + 1). (154) 

The integral in 2 is is trivial at A = and a little complicate at A 7^ but 
still expressible in terms of logarithms of log(l — x); the integral in x in non- 
trivial. Nevertheless one can explicitly check that it is finite and in general 
can be expressed in terms of polylogarithmic functions. In the A — > limit 
and in the on-shell regime p 2 — > — m 2 the full expression strongly simplifies 
and we obtain the expression involving logarithms given in the text. 

With a simple generalization of this method one can compute Feynman 
integrals where the double pole appears, without encountering any problem. 

Moreover, the extension of this method to the integrals appearing in the 
gluon self-energy computation is straightforward. 
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B.3 Feynman parameters integrals 

Finally we report some useful trick to perform integrals on Feynman param- 
eters in particular limits. Suppose one has to compute the integral of some 
function f(x, A) which in the A — > limit is dominate from the x ~ region; 
this is the case for instance for 

f^-W^ir (155) 

where P u (x) is a polynomial a degree v in x, with v < n — 1. This is the 
typical expression we encounter in the computation of self-energy diagrams; 
then one can use the identity 

dxf(x,A)-J dx f(x, A) (156) 

and neglect the second contribution which is subleading at A — > 0. Then the 
first integral can be performed by using the formula 

x m _ B(m + l,n-m-l) /A 2 \ m+1 

o (a 2 + P 2 x) n ~ a 2 ™ [yj ' (157) 

which holds for m > — 1 and n > m + 1. In the case in which the integral is 
dominate by the x ~ 1 region it is sufficient to change the variable x' = 1 — x 
and use the same trick. 



References 

[1] K.G. Wilson and I.G. Kogut, Phys. Rep. 12 (1974) 75 

[2] C. Bagnuls, C. Bervillier, "Exact Renormalization Group Equations. An 
Introductory Review" , |hep-th/ 0002034] 



[3] C. Wetterich, Phys. Lett. 301B (1993) 90 

[4] M.Bonini, M. D'Attanasio, G. Marchesini, Nucl. Phys. B409 (1993) 441 
[5] U. Ellwanger, Phys. Lett. 335B (1994) 364 



46 



[6] C. Becchi, in Elementary Particles, Field Theory and Statistical Me- 
chanics, Eds. M. Bonini, G. Marchesini and E. Onofri, Parma Univer- 
sity 1993, |hep-th/9607T88l ; M.Bonini, M. D'Attanasio, G. Marchesini, 
Phys. Lett. B346 (1995) 87 ; M. D'Attanasio and T.R. Morris, Phys. 
Lett. 378B (1996) 213 



[T 



[9 

[10 

[11 
[12 

[13 
[14 
[15 

[16 
[17 
[18 
[19 



M.Bonini, M. D'Attanasio, G. Marchesini, Nucl. Phys. B444 (1995) 602 

M. Simionato, "Gauge Consistent Wilson Renormalization Group I: 
Abelian case" , |hq>th/9809004| , Int. J. Mod. Phys. A, Vol. 15, p. 2121 
(2000) 

M. Simionato, "Gauge Consistent Wilson Renormalization Group II: 
Non-Abelian case", |hep-th/9810TT7l , Int. J. Mod. Phys. A, Vol. 15, p. 
2153 (2000) 

N. Nakanishi, Prog. Theor. Phys. 67 (1982) 965 

G. Nardelli, R. Soldati, Int. J. Mod. Phys. A5 (1990) 3171 

U. Ellwanger, M. Hirsch and A. Weber, Z. Phys. 69 (1996) 687, Eur. 
Phys.J. CI (1998) 563 

S. Caracciolo, G. Curci, P. Menotti, Phys. Lett. 113B (1982) 311 

J. P. Leroy, J. Micheli and G. C. Rossi, Nucl. Phys. B232, 511 (1984) 

A. Bassetto, G. Nardelli, R. Soldati, "Yang-Mills Theories in Algebraic 
Non-Covariant Gauges" (World Scientific, Singapore 1991) 

D. F. Litim, J. M. Pawlowski, Phys. Lett. 435B (1998) 181 

K. Geiger, Phys. Rev. D60 (1999) 034012 

A. Bassetto, R. Soldati, Nucl. Phys. B276 (1986) 517 

A. Bassetto, I. A. Korchemskaya, G. P. Korchemsky, G. Nardelli, Nucl. 
Phys. B408 (1993) 62 



[20] G. Curci and R. Ferrari, Nuovo Cim. 32 (1976) 151; A. Blasi, N. Mag- 
giore, Mod. Phys. Lett. All, 1665-1674, 1996 



47 



T. Kugo, I. Ojima, "Covariant operator formalism of gauge theories", 
Suppl. of the Progress of Theoretical Physics 66 (1979) 

J. I. Kapusta, Finite Temperature Field Theory, Cambridge Monographs 
on Mathematical Physics, Cambridge Univ. Press, 1989. 

W. Kummer, Phys.Rev. D37, 454 (1988) 

O. Piguet and A. Rouet, Phys. Rep. 76 (1981) 1; O. Piguet and S. 
P. Sorella, "Algebraic renormalization" , Lecture Notes in Physics, Vol. 
m28, Springer Verlag, Berlin, Heidelberg, 1995. 

G. Leibbrandt, "Introduction to Non-Covariant Gauges", Rev. Mod. 
Phys. 59 (1987) 1067; G. Leibbrandt, "Non-Covariant Gauges" (World 
Scientific, Singapore 1994) 

H. Cheng and E. C. Tsai, Phys. Rev. D 36 (1987) 3196 

S. D. Joglekar, A. Misra, Mod. Phys. Lett. A14: 2083-2092, 1999 and 
|hep-th/9912020 



Y. Nakawaki, G. McCartor, "Indispensability of ghost fields and Ex- 
tended Hamiltonian formalism in Axial Gauge Quantization of Gauge 
Fields", lhiFth/0004140 



A. Panza, R. Soldati, "Infrared Singularities in the Renormalization 



Group Flow of Yang-Mills Theories in Axial Gauge" , |hep-th/0006170 



48 



